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Thermodynamics and transport properties of resistance-shunted Josephson junctions are stud- 
ied theoretically in the tight-binding limit Ec/Ej <C 1, where Ec and Ej are a charging energy 
and a Josephson coupling energy respectively. Based on a phenomenological harmonic-oscillator 
model, weak coupling region K = Rq/R <C 1 is analytically studied, where R and Rq = h/(2e) 2 
are a shunted resistance and the quantum resistance. In addition to the effective bandwidth 
2hA e ff, we find that this multi-level system genuinely has a novel crossover at lower energy 
KhA c f{ below which the density of states becomes strongly degenerate. These two energy scales 
control the linear DC responses, optical responses, and nonlinear I-V characteristics. The lower 
energy crossover indicates the existence of a new class of strongly-correlated phenomena beyond 
the framework of the Kondo problem. 
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§1. Introduction 

Small resistance-shunted Josephson junctions are sim- 
plest systems for the purpose of clarifying various macro- 
scopic quantum effects such as quantum t^pjicJjng, quan- 
tum coherence, Coulomb blockade, etcBBBlf Though 
quantum properties of Josephson junctions have been 
studied in many literatures ,a thermodynamics and I-V 
characteristics have not been clarified except for some 
special limits. In this paper, we study specific heat and 
optical linear responses in the weak damping region an- 
alytically. We also conjecture the nonlinear I-V charac- 
teristics from the available results. 

Classical dynamics of the Josephson junctions has 
been studied based on the equivalent circuits shown in 
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Fig. 1. The equivalent circuits for a small Josephson junction 



Fig. 1. The equation of motion is described as 

Io sm<j) = I, 
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CV+ R 



(1) 



where C, R, Iq, I are a capacitance, a shunted-resistance, 
a Josephson critical current, and an external current re- 
spectively. The dynamics of the phase difference (f> gener- 
ates a voltage V across the junction through the Joseph- 
son relation 



(2) 



where $o = h/2e is a unit flux. From (|l|) and (^|), the 
equation of motion is written as 



1 . dU 
Cq+ R q + ^ 



U(q) = -Iq - Ej cos 27T-2- 



(3) 
(4) 



where q = <I>o0/27r and Ej = QqIq/2it. This equation 
can be interpreted as an equation of motion of a dissi- 
pative particle with a coordinate q in the washboard-like 
potential U(q). The mobility fi of this dissipative par- 
ticle is related to the I-V characteristics of the junction 
as 



V = fil. 



(5) 



The Josephson junction has two important energy 
scales: one is the Josephson coupling energy Ej, and 
the other is the charging energy Eq — e 2 /2C. The ratio 
Ej/Ec controls the quantum effects in Josephson junc- 
tions. When the area of the junction S becomes small, 



the charging energy Ec(oc 1/5*) cannot be neglected 
compared with Ej(pc S). Then, the system cannot be 
described by the classical equation. To study dissipa- 
tive quantum systems, a pheiuomenological harmonic- 
oscillator model is introducedfrcP 



H 




Hq = ^c + u(q) ' 



(7) 



where Q is a pseudo-momentum conjugate to q, which 
corresponds to a charge of the junction. The harmonic 
oscillators represent an environment which generates the 
damping term. When the spectral function is taken as 
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the equation of motion (3) is reproduced in the classical 
limit. Here, lo c is a cut-off frequency. 

The ground state of this model has been discussed 
by using perturbative renormalization analysis and dual- 
ity mapping.EP The dimensionless damping strength de- 
fined by K — Rq I R controls the ground-state properties, 
where Rq = h/4e 2 (= 6.5[kf2]) is the quantum resistance. 
For K < 1, the particle is extended in the real space. In 
this case, superconductivity of the Josephson junction is 
destroyed by large phase fluctuations, and the junction 
behaves as an insulator. For K > 1, the particle is local- 
ized in a well, and the superconductivity of the Joseph- 
son junction is retained. Thus, it is predicted that there 
exists the superconductor-insulator transition at K = 1. 
These results have, been supported on a qualitative level 
by experiments.B'Q' 

The I-V characteristics of the resistanre.-shiinterl. htfp. r 
tion have been studied in the literatures.ilWDBB' BBE. 
Generally, it is difficult to study the dissipative model 
for arbitrary ratio Eq/Ej. In this paper, we restrict 
ourselves to the tight-binding limit Ec/Ej <C 1, where 
the system Hamiltonian Hq is well described by a one- 
dimensional tight-binding model 
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in the absence of the external current. Here, A is a hop- 
ping amplitude, and q (c ; ) is an annihilation (creation) 
operator at the site I. Even in this limit, the properties of 
the junction in the presence of the dissipation have been 
studied only for special regions in the whole parameter 
space. So far, there are two regions which allows ana- 
lytical calculations: one is the high-temperature or/and 
strong-damping limit, Elf and the other is the K = 1/2 
solvable lineJlif 



Recently, a fopnulation for K <C 1 has been proposed 
by the authorsJiS Based on this formulation, we can 
calculate thermodynamics and linear responses at any 
temperature. This calculation reveals the independence 
of the low-temperature properties, which have not been 
studied so far. The calculation for the weak damping 
region reveals a novel feature. In the absence of dissi- 
pation, only the bandwidth 2TiA determines the char- 
acteristic energy scale of the junction. Additionally, a 
novel characteristic energy scale KhA appears in the 
thermodynamics and transport properties. This feature 
can be observed most clearly in the weak damping region 
K <C 1, and cannot be observed clearly in the K = 1/2 
solvable case. 

In Ref. |l^, the authors have treated the dissipative 
tight-binding model as an effective model of a dissipa- 
tive carrier in solids, and have not studied applications 
to Josephson junctions. In this paper, by using the re- 
sults of Ref. [TJI and other available references, the aspects 
of thermodynamics and I-V characteristics of Josephson 
junctions are presented in detail. This paper has two 
purposes: (1) we report comprehensive studies on the 
existence of the two different energy scales (especially 
for experimentalists), and (2) we analyze and conjecture 
the junction properties for general values of K beyond 
the solvable cases. 

This paper is organized as follows. The formulation 
of the weak damping region is summarized briefly in 
Sec. |^. The thermodynamics and I-V characteristics are 
discussed in Sec. [| In Sec. ^, the relation to the Kondo 
problem is discussed, and finally in Sec. ^| the results are 
summarized. 

§2. Formulation 

The partition function is formulated from (6) and (|l 
in terms of the path integrals asH 
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where (3 = l/k^T, £j = ±1, and the prime in {£/}' de- 
notes the summation in accordance with the constraint 
6 = 0- The function </>(r) is defined by 

-4 I dw-yW coth(ft/?w/2) 



cosh [w(h(l/2 



smh{hPuj/2) 



(12) 



The expression ( tyj ) may be interpreted as a statistical 
model of classical interacting particles with the potential 
<P(t). 

The weak coupling theory can be constructed by ap- 
plying the Debye-Hiickel approximation to this classical 
partition function.E3' It has not been proved that this 



2 



approximation is valid for the present model. It, how- 
ever, can be shown that this approximation reproduces 
correct results at least in the high-temperature region. 
It also reproduces the correct results in the continuum 
limit (K — > 0, A — ► oo, KA — const.). Hence, we 
expect that this approximation gives correct qualitative 
features. The validity of the Debye-Hiickel approxima- 
tion is discussed in Appendix. 

By the Debye-HiickftLapproximation, the partition 
function is calculated aaU 



Z = 



U(n) 



(13) 



U(n) =n + logT(Kn + 1) - KmJj(Kn + 1), (14) 



Q(n) = K(ij(Kn + l)+j) ) 



2/3HA 



cos a — ne 



-Q(n) 



(15) 
(16) 



where 7 is the Euler constant, T(x) is the gamma func- 
tion, and ip(x) is the polygamma function. The specific 
heat is directly calculated from the partition function. 
The density of states D(tu) is also calculated from Z(0) 
as 



Z(j3) = / duD{u>)e 



-hf)u 



(17) 



where Tiujq is the ground state energy. 

We formulate the optical mobility fi(ui) by the Kubo 
formula. The exact formal expression for fj,(u) is given 
bye 

<J> 2 f°° 

p(u) = -^-Jm dtA{t)e iuJt , (18) 
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where £i,a,a' = ±1, and the prime in {£z,cr, a'}' de- 
notes the summation in accordance with the constraint 
£i + o" + c' = 0. Here, the current-current correlation 
function A(r) can be interpreted as a classical partition 
function with fixed charges at 0, r. Hence, the Debye- 
Hiickel approximation is also applicable in the calcula- 
tion of A(t) in the weak damping region K <C 1. As 
a result, the real-time correlation function A(t) is calcu- 



lated as 
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cos26> smR(t;6), (22) 



where R(t; 9) and S(t; 9) are defined as 

R{t;9) =-KKe-iW\ (23) 



S(t; 
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(24) 



^fc2-(n/3 7 (0)/2^) 2 ' 

7 (0) = 2irKn(9)/hf3. (25) 

Thus, the linear mobility can be calculated at any fre- 
quencies. 

§3. Results 

In this section, we show results mainly for the weak 
damping region K -C 1 based on the formulation of 
Sec. |2[ In this paper, we stress the fact that the two dif- 
ferent energy scales, KhA c ff and hA c g control the junc- 
tion properties. We expect that this new viewpoint is 
helpful to understand the junction properties for general 
values of K and Eq / Ej . 

In addition to the weak damping region K -C 1, there 
exist some regions where analytical -result can be ob- 
tained: the K Y*r\/2 exact solutionis and incoherent 
hopping regime.ElliJ From these results, the thermody- 
namics and I-V characteristics are conjectured beyond 
the weak damping region. 

In the following subsections, we discuss the thermody- 
namics (Sec. |3.ip , the linear DC mobility (Sec. |3.2| ), the 
linear optical mob ility (Sec. |3.3| ), an d th e nonlinear DC 
mobility (Sec. 3.4). Finally in Sec. |3.5| , the I-V char- 
acteristics are derived from a simplified model, and is 
compared with our results. 

3. 1 Thermodynamics 

The specific heat of the dissipative particle has been 
calculated at any temperature in Ref. 15. The low- and 
high-temperature form of the specific heat C is calcu- 
lated for K -c 1 as 



C/h 



2irk B T 

hA cS 

k B T 



2K-2 



12KchA cfi ' 



(fc B T» hA cS ), 
(fc B T< ftAcff), 



(26) 

where the effective hopping amplitude is defined by 
A cff = A(A/w c ) K /( 1 -- R "). The constant c = 
(^Ke*)*/^-^ is close to unity for K < 1. 

In this paper, we focus on strong modification of the 
density of states caused by dissipation. The density of 
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Fig. 2. The density of low-energy states D(ui) in the continuum 
limit for the dissipative case (the solid curve) and the dissipa- 
tionless case (the dashed curve). The density of states for the 
dissipative case has a delta function at uj = 0, which is shown 
schematically by the thick line. 



low-energy states is obtained from thedow-temperature 
expansion of the partition function asH 



D(lu) = VKS{u) + 



1 



(27) 



12VKW 

where W = 2cA c g is the effective bandwidth, and the 
energy is shifted so that the ground-state energy becomes 
zero. This shape of D{lu) is quite different from the one 
in the dissipationless case where D(oj) oc oj~ 1 / 2 . 

We can obtain the whole shape of D(uj) analytically in 
the continuum limit: K — > and A — > oo with KA fixed. 
In this limit, the particle can be treated as a dissipative 
Brownian particle with a damping coefficient.^. =-At:KA, 
which is the only energy scale of the system IL§E3 In the 
continuum limit, the integral (13) is determined by the 
contribution around 8 = 0, and can be evaluated by the 
saddle point method. Thus, the partition function is 
calculated as 

logZ = 2K0A + - \og(2irK) 
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From this expression, the density of states D{uS) can be 
obtained numerically by the inverse Laplace transform 
dl7|). The result takes the form 



(29) 



D(uo) = VK5(lj) + D{uo). 



The excited-state part D(u>) is shown in Fig. 2 by the 
solid curve, while the dissipationless case is shown by 
the dashed curve. The delta function at uj = is also 
schematically drawn by the thick line in Fig. 2. 

If the dissipation effect appeared only in the suppres- 
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Fig. 3. The dispersion relation of the junction. Besides the en- 
ergy scale of the band width 2?iA c ff , there appears a new char- 
acteristic energy scale h'y ~ KhA c ff below which the dispersion 
is strongly modified in the gray region. 



sion of the bare tunneling matrix A to A c g, the band 
dispersion would be given as 



E(Q) = 2hA cS {l - cos(7rQ/e)). 



(30) 



This dispersion, however, does not reproduce the present 
result of D(u)). In particular, the delta function in D{ui) 
at u> = indicates that the low-energy states are degen- 
erate around the band bottom. The qualitative feature 
of the dispersion is shown in Fig. 3, where the dispersion 
relation of relevant excitation is strongly modified in the 
gray region (E < KhA e g). For the moment, it is not 
yet confirmed whether this dispersion is directly related 
to the single particle dispersion relation because single- 
particle Green's function has not been obtained yet be- 
cause of difficulty in the analytical calculation even for 
K < 1. 

As seen in Fig. 3, the dissipative system of K -C 1 has 
two different energy scales: the bandwidth hA c g and the 
damping energy ~ KhA c g. These two energy scales 
appear also in the I-V characteristics as shown in the fol- 
lowing subsections. Note that the existence of these two 
energy scales is not an artifact of the Debye-Huckel ap- 
proximation. Actually, the nontrivial energy scale ap- 
pears in the continuum limit without the Debye-Hiickel 
approximation by mapping the model to the Brownian 
motion problem, and the partition function (28) can be 
reproduced except for the divergent part log K. (See also 
Appendix.) 

As seen in (^7|), the weight of the delta function in 
D{uj) increases as \K for K <C 1. Hence, we expect that 
the macroscopic degeneracy at uj = persists for general 
values of K, and that the weight of the delta function in 
D{uj) increases as K is enlarged. This expectation gives a 
new aspect on the superconductor-insulator transition at 
K = 1. As if approaches unity, the effective bandwidth 
?iA eff = hA(A/uj c ) K/( - 1 ~ K) is suppressed to be zero,© 
and all the states become degenerate. Our result shows 
that a partial degeneracy remains 'even for K < V. 
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Fig. 4. The temperature-dependences of the DC mobility 
fiuc{T) are shown schematically by the solid curve for if < 1, 
and by the dashed curve for K = 1/2. 



Fig. 5. The optical mobility /j.(uj;T) for K < 1. The zero- 
temperature result is shown by the solid curve, and the high- 
temperature result by the dashed curve. 



This degeneracy in the density of states may be ob- 
served by the time-evolution measurement of the volt- 
age, when the external current is large enough to cause 
the Bloch oscillation. If we consider a simplified model 
as discussed in Sec. 3J5, the flat dispersion generates no 
voltage. Hence, the time-dependence of the voltage may 
show the region of zero voltage for some interval. 

3.2 Linear DC responses 

The low- and high-temperature behavior of the DC 
mobility /zdc(^) can be obtained within the linear ra- 
sponse theory for the weak damping region K <C 1 asE£P 



MDc(r) 



1 / k B T 



8n 2 R 



8 \c?iA c ff 
(k B T « HA cS ), 

%A >. 2-2K 

2nk B T / 

{k B T > 7lA e ff). 



+ G(T 4 ) 



(31) 



The qualitative features of /j,uc(T) for K <C 1 obtained 
by connecting high- and low-temperature behavior are 
shown in Fig. 4 by the solid curve. At high temper- 
atures, the mobility increases as the temperature de- 
creases, while at low temperatures, the mobility sat- 
urates to R. The crossover temperature appears at 
T ~ HA cS /k B . 

We expect that this feature of /j.bc(T) is universal for 
< K < 1. In fact, the feature of (|3l|) is valid for the spe- 
cial case K = 1/2, which allows analytical treatment. i$ 
In this case, /j,vc(T) is calculated exactly as 

MDc(T)= J R^f^(aA eff /2fc B T + l/2) ! (32) 
2«b 2 



/ 1 / k B T 
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(k B T < 7iA eff ), 

7T 2 fiA eff 

4fc B T ' 

(k B T » ?iA cff ), 



(33) 



where ip'(x) is the digamma function. The qualitative 
behavior of hdc(T) is shown in Fig. 4 by the dashed 
curve. This result of the K = 1/2 case is consistent with 
our result ©. 

These behaviors of hdc(T) have been interpreted 
physically as follows© At high temperatures, the phase 
<j) slips by 2ir thermally. This phase slip occurs inco- 
herently with a temperature-dependent tunneling rate 
determined essentially by the incoherent hopping in the 
two-level systems. As the temperature decreases, the 
phase-slip picture becomes inadequate, and the quan- 
tum coherence of the system appears. The wave func- 
tion begins to extend in the (f> space. Because of the large 
phase fluctuation, the Josephson junction behaves as an 
insulator, where Cooper pairs cannot tunnel through the 
junction. In this situation, the whole part of the current 
flows through the shunted resistance R, and the voltage 
follows the ohmic law V = IR, which corresponds to 
Mdc = R- 

3.3 Linear optical responses 

In this subsection, we calculate the optical mobility 
/jl(uj;T) within the linear response theory, which can be 
interpreted as the optical responses of the Josephson 
junctions. 

At zero temperature, the optical mobility can be cal- 
culated analytically for K -C 1 at any frequency, and the 
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details of the results have been given in Ref. [IJ. Here, we 
discuss qualitative features of fx(u>) from the analytical 
forms of the limiting regions: 



Re~ 2K ^ 



»(uj;T = 0) 



2T(2K+1) 

7 2 
7 2 + uj 2 



2-2K 



(uj > UJ*), 
(UJ < UJ*), 

(34) 

where 7 is the Euler constant. The qualitative features of 
/z(w) are shown in Fig. 5 by the solid curve. The damping 
frequency is given as 7 = AttKcAcS. Thus, the charac- 
teristic energy scale in /x(w) appears at ^7 ~ KhA c a, 
which is the same order as the energy scale observed in 
the modification of the dispersion. The crossover fre- 
quency uj* in ([m]) is evaluated as uj* ~ 2 1 / K 7. When 
we take K — ► and A — > 00 with KA fixed, the optical 
responses take exact Lorentzian forms (uj* — ► 00), which 
are obtained also by considering the Brownian motion. t3 
As K increases, it is expected that the crossover fre- 
quency uj* in ( |34| ) decreases, and the deviation from 
Lorentzian form becomes clear. At the same time, the 
effective hopping amplitude A e g is suppressed. Partic- 
ularly, in the limit K — > 1, the amplitude A c g vanishes 
and the characteristic frequency 7 ~ K A e g in (|34|) also 
vanishes. For K > 1, the dissipative particle moves only 
by the incohereriLhopping, and the optical mobility takes 
a different formEr 



n(w;T = 0) = R 



8n 2 K 2 A 2 
T(2K + 1)uj 2 



UJ 
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(35) 



The finite-temperature calculation of n(ui\ T) is dif- 
ficult even for K <C 1. For a qualitative discussion, 
however, it is sufficient to consider at high-temperature 
k B T 3> hA c ff, where the optical mobility is calculated for 
arbitrary values of K by assuming the incoherent hop- 
ping of the dissipative particle. The optical mobility is 
obtained for < K < 1 astf 



fi(u;T) = RAttKA cS 
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(uj > Kk B T/h), 

8tt 2 K 2 T(K) 2 f hA cS \ 
T(2K + 1) \2irk B TJ 
(uj < Kk B T/h). 

The qualitative behavior is shown by the dashed curve 
in Fig. 5. The high-frequency side is temperature- 
independent, and agrees with the zero-temperature re- 
sult in the weak damping region if « 1. This shows 
that the high frequency side is well described by the in- 
coherent hopping of the dissipative particle even at zero 
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Fig. 6. (a) The obtained I-V characteristics for < K < 1, and 
(b) the exact solution at K = 1/2. The zero-temperature results 
are shown by the solid curves, and the high-temperature results 
by the dashed curves. The Ohmic responses V = IR are also 
shown by gray lines for references. 



temperature. The low-frequency side is temperature- 
dependent, and behaves as T 2K ~ 2 . The crossover fre- 
quency is given asw~ Kk B T/Ti. 

When the temperature decreases, the low-frequency 
temperature-dependent part in (36) becomes inadequate 
below a crossover temperature. This crossover temper- 
ature is estimated for K 1 as T ~ ?iA c ff/fcB where 
the characteristic frequency in (36) takes the same order 
(7 = AttKA cS - Kk B T/h) as that seen in (|J). This 
result is consistent with that of the DC mobility, where 
the crossover temperature is given also as T ~ hA e s/k B . 

3.4 Nonlinear responses 

Theoretically, it is a difficult task to study nonlin- 
ear responses beyond the linear response theory, because 
nonequilibrium states must be treated properly. In this 
paper, general properties of the I-V characteristics are 
inferred from the results obtained in some limits. 

At zero temperature, the nonlinear mobility /xj\jl = 
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V/I can be calculated for large I asl 
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where e = 7<I>o is the energy difference between neigh- 
boring wells in the periodic potential U(q). For large /, 
the voltage V = /xnl-^ is proportional to I 21 ^^ 1 . On the 
other hand for small /, the mobility /xnl is expected to be 
R far < K < 1 from the renormalization group anal- 
ysisU By connecting these two limiting cases smoothly, 
the I-V characteristics are derived as in Fig. 6 (a) by the 
solid curve. 

The crossover current is determined by the cross point 
of the results of the two limiting cases. Particularly for 
K -C 1, the crossover current is e ~ KhA e s, while the 
characteristic voltage is V ~ IR ~ ft,A e g/e. Thus, also 
in the I-V characteristics, two different energy scales, 
hA c f{ and KhA^g appear clearly. 

The high-temperature limit k B T ^> hA c g can also be 
calculated analytically for < K < 1 as 
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(40) 



T{2K + \) \27rk B T 
(e< fc B T). 

The schematic J-V characteristics are shown also in 
Fig. 6 (a) by the dashed curve. At large current (e = 
3> k B T), the mobility is temperature- independent, 
and behaves as I 2K ~ 2 , while at small current (e <C 
k B T) it becomes temperature-dependent. Note that in 
this regime the nonlinear DC mobility is obtained from 
the linear optical conductivity as fJ.Nh(£]T) — m( w ~> 
e/h-T). 

We expect that the obtained I-V characteristics are 
valid for any values of K. This conjecture is easily 
checked by considering the exact solution of the K =1/2 
case given byll3 
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The I-V characteristics in the K = 1/2 case is shown in 



Q 



i. 



(42) 



Fig. 6 (b) for k B T <§; TiA c s by the solid curve, and for 
k B T ^> TiA e fi by the dashed curve. This I-V character- 
istics are consistent with the ones obtained in Fig. 6 (a) 
when we take K — 1/2 in our present result. 

3.5 Discussion based on a simplified model 

The above results can be understood intuitively in 
terms of the Coulomb blockade and the Bloch oscillation 
by a simple analysis based on the equation of motion 

Ye 

R 

which is a simplified version of the one analyzed in Ref. || 
in detail. Here, Q and I are the charge in the junction 
and the external current, respectively. The second term 
in (|2|) denotes the normal current through the shunted 
resistance. The voltage is calculated by V x = dE(Q) / dQ, 
where E(Q) is the dispersion relation given by (|30|). 
The kinetic energy of the particle originates from the 
Coulomb energy, and C c s — e 2 /TiA c r plays a role of 
an effective capacitance. The Josephson effects appear 
from the fact that the pseudo-momentum Q is folded 
into the region — e < Q < e. For the tight-binding limit 
Ej/Ec ^ 1, it is allowed to focus only on the lowest 
band. We note that ( p2| ) is justified only when Q is a 
classical variable and quantum fluctuations are negligi- 
ble. 

By solving the equation ([42]), we find that the junction 
properties change at a critical current I* = 2irhA Q g / eR. 
For I < I* , the charge converges to a finite value (Q = 
0). In this regime, a Cooper pair cannot tunnel through 
the junction, because it costs additional Coulomb energy 
for the Cooper-pair to tunnel. This behavior is easily 
understood as the Coulomb blockade phenomena. 

For the opposite case I > I* , the charge Q continues 
to increase until Q becomes e, and then the value of Q 
jumps to — e by tunneling of a Cooper pair. This pro- 
cess iterates permanently. In this regime, there appears 
voltage oscillation withj-the frequency / = I/2e known 
as the Bloch oscillation^ 



From the simplified model (|42j), the voltage averaged 
over the period T = 1/f is given as 

V = J q T ^V x , (43) 

and is calculated as a function of I. The result is shown 
in Fig. 7. For I < I* , the I-V characteristics are ex- 
pressed exactly as V = IR, because no Cooper pair 
can tunnel through the junction by the Coulomb block- 
ade. For I > I*, the junction retains superconductivity 
by opening the channel of the Cooper-pair tunneling, 
and the voltage is reduced. The whole shape of this 
I-V characteristics is similar with the one obtained in 
the previous subsection. Particularly, the characteristic 
voltage V* and current e* = J*$o are consistent with 
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Fig. 7. The I-V characteristics calculated from the simplified 
model Ji^). The current and voltage is normalized by I* = 
2nHA c ff /eR and V = 2nhA c f{/e, respectively. The dashed line 
denotes the Ohmic response V = IR. 



the one obtained in the last subsection (V* ~ frA e g/e, 
e* ~ KHAeff). Thus, the two energy scales hA e g and 
KhA c g appear also in this simplified model. 

We, however, should note that the simplified model 
does not reproduce correct results at a quantitative level 
due to the absence of quantum fluctuations of Q. The dif- 
ferences can be pointed out as follows: (1) the asymptotic 
behavior for I I* oi the simplified model is given by 
V oc 1/1, while the exact calculation gives V oc I 2K ~ 1 , 
(2) the exact solution for K = 1/2 suggests that there is 
no cusp in the I-V characteristics for any K in contrast 
to the result obtained from the simplified model, and (3) 
the exact solution for K = 1/2 also clarify that even for 
small I the superconductivity of the junction is partially 
retained as seen in the deviation of the I-V characteris- 
tics deviate from the V = IR line. Thus, it is clear that 
the simplified model cannot be used for the quantitative 
discussion. 

§4. Relation to the Kondo problem 

In this section, the dissipative tight-binding model is 
discussed in relation to the Kondo problem. Although 
this discussion is not helpful in analyzing unsolved prob- 
lems about the dissipative tight-binding model, it pro- 
vides us with a new perspective, which may contribute 
to a better understanding of the junction properties. 

First, we consider dissipative systems with a double- 
well potential. When the systems are truncated to two- 
level systems, the Hamiltonian is obtained as 



hA 



-E 



2m a 




where a x and a z are the Pauli matrices. This model has 
been studied as the simplest model describing destxi. 
tion of quantum coherence by surrounding medium. □' I 
In this paper, the dissipative 'multi-level' system is 



studied. One may regard this model as an extension of 
the dissipative 'two-level' system. There, however, exists 
a critical difference between these two models. In dissi- 
pative two-level systems, the dissipation is introduced 
to the 'gaped' system with a gap hA, while in multi- 
level systems, to the system with a 'gapless' band with 
a bandwidth 2hA. Though these two models show es- 
sentially the same behavior at high temperature, the 
low-temperature properties are quite different. Actu- 
ally, density of states for dissipative two-level systems 
has no characteristic feature at the energx,scale K%A e g 
in contrast with the multi-level systems.EZP The energy 
scale KhA c g appears in the dissipative two-level systems 
only as a characteristic relaxation time in the correlation 
function and in_^the evolution of the expectation value 
P(t) = (a,(t))B _ 

Even in this dissipative two-level model, it is a nontriv- 
ial task to study the system properties such as thermo- 
dynamics and real-time correlation functions. Actually, 
it has been proven that the low-energy behavior of this 
dissipative two-level model is equivalent to—the Kondo 
model with anisotropic exchange couplings:^ 



SySy). 



(45) 



k.a 



Here, Ck.a (ct CT ) is an annihilation (creation) operator 
of a conduction electron, Si is the local spin operator, 
and Si is the spin operator of the conduction electrons at 
the impurity site. The equivalence between the Kondo 
model and the dissipative two-level system is obtained by 
identifying the coupling constants Jm and J± withH 



A 2 7T/?J| 

— = p J± cos — — 



K = 1 - 



P J \\ 



(46) 
(47) 



where p is the density of states at the Fermi surface. 
Unfortunately, the weak damping region K <C 1 cannot 
be related to the Kondo problem exactly, since the abojse 
relations become valid only in the limit pJn, pJ± <C 1.E3' 
The properties of the dissipative two-level model can 
be understood in terms of the Kondo problem as fol- 
lows.B 

In the dissipative two-level systems, the effective hop- 
ping amplitude A e ff determines the energy scale of the 
system. At high temperatures k^,T 3> hA e g, the coher- 
ence of the system is destroyed by dissipation, and the 
energy levels become degenerate. The particle is local- 
ized at one well, and hops to the other incoherently. At 
low temperatures k^T <C foA c g, the coherence of the 
system is retained, and the energy splitting is given by 

TlAeff. 

In the (anisotropic) Kondo model, the Kondo temper- 
ature Tk determines the energy scale of the system. At 
high temperatures T 3> Tk, the local impurity spin fluc- 
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tuates thermally, and the magnetic susceptibility of the 
local spin \ follows the Curie law. At low temperatures 
T <C Ik, the local impurity spin forms a singlet state 
coupled with conduction electrons, and \ becomes inde- 
pendent of the temperature. 

In the both cases, at a characteristic temperature 
(hA e $ or k^Tx), the behavior of the system changes 
drastically. This common feature actually comes from 
the equivalence of these two models-EScEP 

From this viewpoint, let us come back to the results 
obtained in the previous section. In the dissipative multi- 
level model, there are two different energy scales, which 
are seen clearly in the calculation for K -C 1. One is 
the effective bandwidth 27iA e g, which can be related to 
the Kondo temperature 7k. This energy scale appears 
as a crossover temperature in the DC mobility, and as 
a characteristic voltage V ~ hA e s/2e. The other is the 
new energy scale fvy ~ 7T?iA c ff , which can be seen in the 
density of states D(u>) (see Fig. 3), and in the optical 
mobility fJ>(u>). This energy scale is characteristic of the 
multi-level model, in which the dissipation is introduced 
to a gapless system. This situation cannot be realized in 
the Kondo problem. 

In this sense, the resistance-shunted Josephson junc- 
tion system offers a new type of nontrivial strongly- 
correlated problems. Though it has some similarity, it 
also has an essential difference from the Kondo problem. 
If properly designed, Josephson devices may show such a 
novel phenomenon originating from strongly correlation 
effects. 

§5. Summary 

In this paper, thermodynamics and linear responses for 
the limiting region (K <C I and Ej/Ec ^> 1) have been 
studied analytically for all temperatures. The nonlinear 
I-V characteristics are also discussed from the available 
results. 

The density of low-energy states (cj < KA e g) is 
strongly modified. This modification indicates that the 
dispersion relation of relevant excitations is strongly al- 
tered at the band bottom. The DC mobility (J,t>c(T) — 
V/I behaves as T 2K ~ 2 at high temperatures, while it 
saturates to R below a crossover temperature T ~ 
hA e s/kte. At zero temperature, the optical responses 
follow a Lorentzian form with a width 7 = 4TtKA e f[ at 
low frequencies, while they behave as lq 2K ~ 2 at high fre- 
quencies. Above a crossover temperature T ~ 7iA c ff//cB, 
the optical mobility becomes temperature-dependent on 
the low-frequency side. The nonlinear I-V character- 
istics for < K < 1 change at a crossover current 
I* ~ KhA c ff/$Q. For I < I*, the system behaves like 
an insulator following the ohmic law V = IR, while for 
I > I* the junction retains the superconductivity, and 
behaves as V ~ P K ~ l . These features are reproduced 
qualitatively by a simplified model. Among others, this 



system shows two different energy scales. One is the ef- 
fective bandwidth HA c ff, and the other is KhA c s. The 
latter energy scale appears as the consequence of strong 
coupling to the dissipation and is seen in the density 
of states, optical mobility and I-V characteristics. The 
Josephson devices can be recognized as novel strongly- 
correlated systems. 

In this paper, the junction properties have been stud- 
ied mainly in the weak damping region K <C I and the 
tight-binding limit Ej/Eq ^> 1. These two conditions 
can be realized easily in experimental situations. In the 
tight-binding region, however, the effective bandwidth 
SA e ff must be small (typically hA e g/kB < I [mK]). 
Then it becomes difficult to lower the temperature be- 
low hAes/ks- For experimental feasibility, the region 
Ej ~ Eq would be more fruitful, where the qualitatively 
similar feature predicted in this paper can be observed. 
Actually, also in the opposite limit Ej/Eq -c I, two 
different energy scales appear: the effective band width 
A = e 2 /2C determines the characteristic temperature of 
the linear mobility, while the damping energy 7 ~ KA 
determines the characteristic current e = 7$o of the non- 
linear mobility.™ The analyses in this paper may con- 
tribute to the understanding of quantum effects in the 
Josephson junctions also in those parameter regions. 
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Appendix: Validity of the Debye-Hiickel Ap- 
proximation 

The concept of this approximation has first been con- 
sidered by Debye and HuckeLKiP Later, Mayer has shown 
that the Debye-Hiickel theory can be reproduced from 
the ring approximation in the classical cluster expan- 
sionO Since then, the corrections to the Debye-Hiickel 
approximation have been studied for the problem of in- 
teracting plasma. 

Let us consider the 3D classical Coulomb gas with the 
particle density n. In the Debye-Hiickel approximation, 
the two-body distribution function is calculated as 



Here, T is the temperature, and 4>T>{r) is the screened 
potential 

Mr)=e 2 ^—, (A-2) 
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where 1/kd is the Debye screening length given by 



T 



2\ 1/2 



(A-3) 



Note that the two-body distribution function can be 
written in a simple form as 



P2 (r) = n 2 exp 



(A-4) 



where x — K^r is the dimensionless spatial coordinate, 
and £ = e 2 KT>/T is the plasma parameter. Beyond the 
ring approximation, it is possible to sum-the-graphs up to 
the order of e 2 in the cluster expansion Bj'EfP The result 
is given as 



p 2 (r) = n exp — e e fix) 

x 



(A-5) 



where f(x) is a complex function of x. From this result, if 
the series of e is convergent, then the ring approximation 
is justified for e ^J.. Also the results b y the Monte 
Carlo simulations^' Et have supported Eq. ( A-5 ). Thus, 
the Debye-Hiickel approximation is thought to be valid 
in the weak coupling limit e <C 1 for the 3D classical 
Coulomb gas. 

One may suspect that the Debye-Hckel theory is not 
applicable to one-dimensional systems. Here, let us 
discuss the ID classical Coulomb gas described by the 
Hamiltonian 



H = - p? - 27re2 Y °i°Mi - <ij i 



i<j 



(A-6) 



where pj, qj, and Oj — ±1 are the momentum, posi- 
tion, and sign of the j-th charged particle. Among 2N 
particles, N particles are charged positively and N neg- 
atively. Fortunately, this model can be solved exactly.EZP 
We take the dimensionless parameter as 7 = P/2-7re 2 , 
where P is the pressure. In this ID model, the system 
behaves as ideal gas in the high-pressure limit 7 — > 00 
on the contrary to the 3D case. The energy per particle 
u is calculated for 7 3> 1 in the expansion form as 



T 



^1 

2 2 



1 



1 



1 



(27)V2 87 128V27 3 / 2 



(A-7) 

It can easily be checked that the Debye-Hiickel approx- 



imation gives the exact result (A-7) up to the order of 
I/7 1 / 2 . Thus, this approximation is also applicable to 
the ID classical Coulomb gas for the weak coupling re- 
gion 7 > 1. 

From the above discussions, the Debye-Hiickel approx- 
imation is thought to be applicable to general classical 
plasma models in the condition 



1/«d > I, 



(A- 



where I is the averaged particle distance. This condi- 
tion is in fact reasonable from the physical viewpoint 



because the Debye-Hiickel approximation, an effective 
field approach, is expected to become reliable when many 
particles are interacting each other within the screening 
radius. To clarify further the validity of the ring approx- 
imation, other methods such as Monte Carlo simulations 
would be helpful. 

To study the dissipative tight-binding model, ID clas- 
sical particles interacting with the potential 4>{t) in the 
r direction have been considered. The potential is for- 
mulated as 



hp 



1W„ 



2itK 



Mm 



(A-10) 



for < \bj. 
uj m > uj c . The condition (A 



whi le t he potential vanishes for 
is evaluated for the dis- 
sipative tight-binding model as follows. In the Debye- 
Hiickel approximation, the screened potential ip(iuj m ) is 
calculated as 







1 - ruf)(iu) m )/h(3' 



(ATI) 



The screening effects appears when the second term of 
the denominator in ( |A-11 ) is dominant. Hence, the De- 
bye screening length 1/kd is determined from ri(f>{u> m — 
Ko)/Tif3 ~ 1, and the condition ( |A-8| ) corresponds to 
the weak-damping condition K <C 1. Thus, we expect 
the ring approximation gives reliable results in the weak- 
damping region. 

It can be proven that the Debye-Hiickel theory repro- 
duces correct results in_the high temperature and/or 
strong damping region.H The Debye-Hiickel approxi- 
mation also reproduces the zero-temperature mobility, 
wjj|ij|i is obtained by the renormalization group analysis 



Hdc(T = 0) = R. 



(A12) 



Further, in the continuum limit (K — > 0, A — * 00 with 
7 = 4irKA fixed), it can be proven that the Debye- 
Hiickel approximation gives corrects results for the op- 
tical mobility^ Also for thermodynamics, it can be 
checked that the energy calculated in the Debye-Hiickel 
approximation 



E — const. H 



2^ 



1 



2tt 



np 1 \ 

(A-13) 

agrees with the energy of a dissipative Brownian particle 
obtained in the continuum limit as 



h 

2^ 



duj- 



uj + 7- 



coth(wS,/3/2). 



(A-14) 



These provide further supports for the reliability of the 
Debye-Hiickel theory in this problem. 
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